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Abstract 

Quasinormal modes of asymptotically AdS black holes can be interpreted as poles of retarded correlators 
in the dual gauge theory. To determine the response of the system to small external perturbations it is 
not enough to know the location of the poles: one also needs to know the residues. We compute them 
for R-charge currents and find that they are complex except for the hydrodynamic mode, whose residue is 
purely imaginary. For different quasinormal modes the residue grows with momentum q, whereas for the 
hydrodynamic mode it behaves as a damped oscillation with distinct zeroes at finite q. Similar to collective 
excitations at weak coupling the hydrodynamic mode decouples at short wavelengths. Knowledge of the 
residues allows as well to define the time scale Th from when on the system enters the hydrodynamic 
regime, restricting the validity of hydrodynamic simulations to times ? > Th- 
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I. INTRODUCTION 



In recent years a new paradigm concerning the high temperature behaviour of QCD has been 
established: the strongy coupled Quark-Gluon Plasma (sQGP). Experimental results from heavy 
ion collisions at RHIC indicate that QCD at temperatures around 2Tc is strongly interacting, in 
spite of being in a deconfined phase, and thus rendering perturbative computations not suitable for 
describing it. While static properties of strongly coupled gauge theories at finite temperature can 
be readily analyzed on the lattice, the study of out-of-equilibrium phenomena faces considerable 
difficulties. In the last years the AdS/CFT correspondence has emerged as a useful tool to un- 
derstand analytically the dynamics of non-Abelian gauge theories in the strongly coupled plasma 
phase. 

Perturbative computations in gauge theories show that, at finite temperature, zero momentum 
correlators of gauge-invariant operators have an inifinite set of discrete poles in the complex fre- 
quency plane These are responsible of the exponential decay of correlators with time, so 
describe the dissipation of small perturbations in the plasma. According to the AdS/CFT conjec- 
ture, the asymptotically AdS black hole is dual to the plasma phase of the strongly coupled 3\f = 4 
gauge theory [13D. Correlation functions in the thermal gauge theory can be computed using classi- 
cal solutions of fields living in the gravity dual \^. If a field is excited in the presence of a black 
hole, the energy of the fluctuation will be lost inside the horizon and eventually the final state will 
be a larger black hole with no fluctuations. This process is described by the quasinormal spectrum, 
that was first computed in [|5|] for black holes in asymptotically Anti de Sitter (AdS) spacetimes. 
Quasinormal modes (at fixed momentum) exist only for a discrete set of complex frequencies that 
can be identified with the poles of the retarded Green's functions in the dual field theory fla lZO . 
Therefore, quasinormal modes describe dissipation processes in the plasma. 

Hydrodynamic modes like diffusion, shear or sound modes are also described by quasinormal 
modes in the gravity dual. One of the most interesting results has been the derivation of a universal 
bound for the shear viscosity to entropy ratio ^ > 4^ [8]. It has also been argued that this 
value is relevant to the description of heavy-ion phenomenology at RHIC [9] . One way to derive 
it is to compute the lowest quasinormal mode in the retarded two-point correlator of the stress 
tensor. It is an example of a hydrodynamic mode; a special case of quasinormal mode whose 
frequency vanishes in the zero momentum limit [llOn . Another interesting application used the 
spectral function of R-charge current correlators to derive photon and dilepton production rates in 
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the strongly coupled plasma fll^. 

It is helpful to remember the situation at weak coupling for correlators of field components: 
the poles of the retarded Green's functions correspond to quasiparticle excitations UM- A hard 
scale of order T can be distinguished from a soft scale of order gT. At the hard scale the relevant 
excitations are the hard partons, i.e. the quarks and gluons. At the soft scale however there are 
quasiparticles corresponding to the dressed partons and collective excitations, e.g. the longitudinal 
plasmon mode. Both types of poles can be distinguished by the behaviour of their residues at short 
wavelengths: the residues of the particle poles scale like q^^ with the momentum q whereas the 
collective excitations show exponentially decaying residues of the form exp{ — aq^/g^T^) with a 
mode dependent constant a [12]. To gain a better understanding of the quasinormal frequencies 
that appear in the holographic model of the strongly coupled plasma it is of utmost importance to 
study the behaviour of their residues. In particular we will concentrate on the transverse and longi- 
tudinal R-current correlators. The corresponding quasinormal frequencies and spectral functions 
have been studied before in [[ibI] and [ 14,ll5], respectively. 



Knowledge of the residues at these frequencies is necessary to compute the response of the 
system to an external perturbation. It also allows to define the time Th from which on a hydro- 
dynamic approximation is valid. Since hydrodynamic simulations of the evolution of the Quark- 
Gluon plasma are at the heart of the new sQGP paradigm, it is of utmost importance to know the 
hydrodynamic time scale Th, which in turn requires knowledge of the residues of the quasinormal 
modes. 



n. LINEAR RESPONSE 

The response to an external perturbation represented by the source x) is given by 

(4>(?,X)) = -J dTd3^GR(?-T,X-^)7(T,^) , (1) 

where Gr is the retarded two-point function. This expression is valid when the energy of the 
perturbation is negligible compared to the total energy of the system. If we choose a perturbation 
that is localized in time, i.e. we strike the medium once at time t = and with a periodic profile 
such that j{t,x) = 5(f)cos(qx), we find 

(4>(f,x)) = /0(Ocos(qx)£7?„(q)e-'"«('i>-r„(q> ^ 

n 
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where we have assumed that the only singularities of the retarded Green function on the complex 
frequency plane are poles at the values = t2„ — zT„ and 7?„ are the residues. Stability demands 
that all poles lie in the lower half plane. The contribution of a particular pole is proportional to the 
residue; hence their importance. 

There is a restriction on the location of the poles and on the values of the residues. This 
restriction can be derived either from the condition that the response ^ has to be real or from 
symmetry under time reversal p(a),q) = — p(— (0,q) where p((D,q) = — 2ImGR(ft), q). Both 
conditions lead to the result that the poles come in pairs and that the residues at a pair are related 
according to 

= Rn = -K. (3) 

There can also be unpaired poles lying on the imaginary axes and in this case their residues have 
to be purely imaginary. 

We are able to compute the location of the poles and the value of the residues using the gravity 
dual. They are determined by the values of the frequency where quasinormal modes exist, and 
by the shape of these classical solutions. In our numerical calculations we have checked that the 
stability and time reversal relations indeed hold for the quasinormal modes. Analytic proofs can 



be found in [|l7l] 



Of particular interest are the hydrodynamic modes, whose frequencies vanish as the momentum 
goes to zero. At sufficiently small momentum the hydrodynamic mode has the smallest imaginary 
part of all poles and therefore dominates the long time behaviour. Knowledge of the residues 
for the quasinormal modes allows us to define a hydrodynamic time scale Th» from which on the 
hydrodynamic mode dominates the response|fl^. From Q we find 

l0g|7?H|-l0g|i?l| 

= :^ :^ • (4) 

i H-i 1 

Only for times larger than Th we can expect hydrodynamics to be a good approximation. For 
shorter times more and more higher poles will contribute. 



ni. AN EXACTLY SOLVABLE EXAMPLE 



In general the quasinormal modes of the five dimensional AdS black hole can be determined 
only numerically. There are however some special cases where the wave equations simplify and 
can be solved analytically. Two such cases are the wave equations for the gauge invariant variables 
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corresponding to longitudinal and transverse vector field perturbations El = Mq + (D q ■ A/^, E'j = 
(oAj. The AdS black hole background is given by 



d^' 



^(-/Wdr2 + dx^) + f|J 



(5) 



with /(r) = {I — Tq/t^). The Hawking temperature Th = ro/nL?- is the temperature T of the dual 
field theory. In the following we will use the coordinate x = I — r^/r^ such that the horizon sits at 
X = and the boundary atx= 1 . The equations of motion for these gau ge in variant combinations 
of the vector field perturbations with frequency CO and momentum q are lll3ll 



E'i + ^-E'r + 



/ 



2 f/ 











(6a) 
(6b) 



where we have introduced the dimensionless frequency and momentum 2nT{)xi,q) := ((0, q). Ac- 
cording to the holographic dictionary the vector field in AdS acts as a source for the R-charge 
currents of the SU(4) R-symmetry in the dual gauge theory. Denoting the retarded correlator 
of two currents as G^v one can write for the non-vanishing components jlBIl 

rn^((0,q) 
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Gtt = n^(w,q) 



(0' 
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GtL 



-qCO 



n^((0,q), Gll 



CO" 



rn^(«,q) 



G)^ - q^ (o^ - q" 

The correlators are therefore defined by the transverse and longitudinal polarization tensors 
n^''^((i),q). 

At zero momentum q = longitudinal and transverse components become indistinguishable 
from each other and obey the same wave equation. Moreover at q = 0, equations Q can be solved 
exactly and the polarization tensor was found to be [|l5n 



n(ro) 



in) + to 



¥\ ^V^tt' ]+¥ 



[l+i) 



-Ki 



(8) 
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where i// is the digamma function and the number of colors. The locations of the poles and their 
residues are 



ft), 



± 



2nTn{±l-i) , = '-^N^T'n^{T\-i) 



(9) 



with n eN. The residues fulfil the reality conditions Q. At first glance it looks strange that the 
residues grow quadratically with the mode number. If we want to write the retarded Green's func- 
tion as a sum over the poles this leads to a divergent expression. However, we have to remember 



that the poles determine the Green's function only up to an analytic part and that this analytic part 
can be infinite. In order to obtain a well defined pole representation one has therefore to substract 
a conveniently chosen analytic part. This is made explicit in the well-known pole representation 
of the digamma function 

¥{x) = -Te - £ ( \— - -) , (10) 

where is the Euler-Mascheroni number. Using this we can find a pole representation of the 
polarization tensor as 

U{xv) = ^l-irv-2rErv^ + rv'y-( \ - + \ -]] . (11) 

The linear term in ro can be understood as the zero momentum limit of the diffusion pole. This sug- 
gests that the retarded Green's function at non-zero momentum would be a meromorphic function 
of frequency and momentum and admit an expansion 

where i?„(rD, q) are analytic functions. We will do a numerical evaluation of the residues at the 
pole i?„(rD„(q), q). This implies that our method does not disentangle the frequency from the 
momentum dependence. 

A rather interesting point is to ask what is the spectral line produced by a single isolated quasi- 
normal mode. Due to the fact that the residue is complex the line spectrum of an isolated quasi- 
normal mode does not have the shape of a Lorentzian function. Writing Y = lm{R)/Re{R), we 
find 

, , r-Y{rv-a) 

For non- vanishing imaginary part of the residue the peak is higher than the one of a Lorentz curve 
with y = and it is shifted to lower frequencies if 7 > or to higher frequencies if 7 < 0. In the 
case at hand where the T/Q.^ I the quasiparticle approximation is not valid and the line spectrum 
of an isolated quasinormal mode does not approximate the full spectral function in any range of the 
frequency. In cases where r/Q.<^l and Y different from zero it would show up as an asymmetry in 
the form of the observed resonance in the spectral function and a correct interpretation is possible 



only if the residue is known [|20|l 



The expression found as a sum over poles and the exact one differ only by some contact terms, 
coming from the analytic pieces we have neglected. Neglecting the contact terms, the Green's 



function diverges at ? = as t~^, coinciding with the zero-temperature result. This shows that the 
sum over poles is enough to recover most of the dynamical information and only the contact terms 
require a more detailed analysis. In principle they should be related to the Schwinger terms of the 
r = theory. 



IV. NUMERICAL COMPUTATION 



According to the holographic dictionary the retarded Green's function can be computed as 
the ratio of the connection coefficients that relate the local solution at the horizon with ingoing 
boundary conditions to the non-normalizable (^1) and normalizable (B) solutions at the boundary. 
Defining (a) := (r,L) as the two components 



^(a)W=-^(«)^(a)W+S(a)^(a)W- 



(14) 



The retarded polarization tensors are given then by [|l3h 

A^2 j2 3 



n(«) 



8 A 



(«) 



The quasinormal modes are normalizable solutions where the connection coefficient A 



(a) 



(15) 



0. 



The polarization tensor follows from demanding that the solution is smooth at a matching point in 
the interior of the interval (0, 1), 



A 



(a) 



^ {a)' (a) 



(16) 



We have computed the Frobenius series up to order 50. Matching the series expansions, we see 
that the ratio (fT6l) remains constant for a fair interval in the radial coordinate. We have chosen 
X = Q.53 to evaluate the ratio and have checked that the spectral function agrees with previous 



numerical (for non-zero momentum) and exact (for zero momentum) results yj 
can be computed as 



ISfl. The residue 



R 



(a) 
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3m 




ro=ro„- 



n -1 



(17) 



Our results for the residues (normalized by [{nN^T^n'^) /4]^^) of transverse and longitudinal vec- 
tor fluctuations are plotted in Figure [IJ The longitudinal fluctuations show interesting behaviour 
related to the diffusion mode. The peaks and dips in Figure \T\ appear roughly at the locations 
where the hydrodynamic mode [12 ill crosses the imaginary part of the quasinormal frequency. The 
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hydrodynamic mode behaves in a completely different way. The diffusion pole quickly moves 
towards negative imaginary frequencies, while the residue first grows according to hydrodynamics 
and later goes over into a damped oscillation. Numerically the zeroes coincide with ro = —in. 
On the other hand, the location of the other quasinormal modes remains fairly constant until the 
diffusion pole reaches an integer value that coincides with the imaginar y part of the quasinormal 
pole at zero momentum. After this point, the quasinormal pole starts approaching the real axis. 
This first happens at = 1/2, when the diffusion pole is at ro = —i. For larger momentum, the 
long time behaviour of the fluctuations will be dominated by the first quasinormal mode and not 
by the hydrodynamic mode. We also find some interesting analytic structure related to the zeroes 
of the residue of the hydrodynamic pole (see Figure [2l). When ro = —in, the exponents of ingo- 
ing and outgoing solutions at the horizon differ by an integer. The ingoing solution should then 
develop a logarithmic term since it has the lower exponent. However, for very special values of 
the parameters, i.e. the momentum q, it can happen that the coefficient of the logarithm vanishes. 
We find numerically that the zeroes of the residues of the hydrodynamic mode coincide with these 
special values at (tt),q2) = (-/, 1/2), {-2i,V3-l), {-3i,V6-3/2) and (-4/, 1.4436). 



V. DISCUSSION 



The analytic structure found for the R-charge correlator at finite temperature is very interest- 
ing. At zero momentum, an infinite set of evenly spaced poles on the complex frequency plane is a 
generic feature of correlators. At weak coupling, they lie on the imaginary axis M\, but according 
to the AdS/CFT computation they move away at strong coupling. The right UV behavior, that 
is analytic at zero temperature, is recovered after summing over all the poles. The value of the 
residues is crucial, and especially the vanishing of the hydrodynamic mode at zero momentum. 
When we consider momentum dependence, the analytic form of the diffusion residue in the hy- 
drodynamic approximation ~ will induce a singularity in the correlator at short times, on top 
of the usual UV singularity. This clearly cannot be the right answer, and we can understand this as 
a limit in the validity of the approximation. The damped and oscillatory behavior of the hydrody- 
namic residue can cure this problem, making the contribution from the diffusion mode smooth at 
short times. The hydrodynamic mode effectively decouples for momenta q> I and in this respect 
behaves as the collective modes present at weak coupling. In this sense, our results go beyond the 
hydrodynamic approximation. In principle, we expect that other hydrodynamic modes appearing 
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Figure 1 : (Top) Real and imaginary parts of the residues for the first four quasinormal modes in the trans- 
verse component Ej (bottom) Idem for longitudinal component Ei^. The «^ scaling is necessary to recover 
the asymptotic behaviour of the spectral function at large frequencies. Close to the crossing with the dif- 
fusion mode q ~ 1, the residues of the longitudinal component present peaks. The residues grow with 
momentum, this is also reflected in the growth of the spectral function. 

in the stress-tensor two-point functions, the shear and the sound mode, to have a similar collective 
mode behavior. 

We observe that as the momentum increases, diffusion becomes less important and other col- 
lective excitations of longitudinal modes describe charge density fluctuations. In contrast with 
the dressed partons of the weak coupling regime, they do not decouple at high momentum. The 
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Figure 2: In the lower figure, we show the real and imaginary parts of the quasinormal frequencies in the 
longitudinal channel and the value of the frequency for the diffusion mode. The quasinormal frequencies 
remain fairly stable as momentum increases, until there is a 'crossing' with the diffusion mode (it reaches a 
special value ro = —in). Then, there is a qualitative change in the behaviour of the quasinormal frequency 
that starts approaching the real axis. This behaviour can also be oberved in previous computations of quasi- 
normal frequencies 12]. The residue of the diffusion mode, shown in the upper figure, behaves according 
to hydrodynamics ~ for small momentum. For larger momentum, it shows an oscillatory and decaying 
behaviour. The zeroes of the residue coincide with the 'crossing' values. Quasinormal frequencies in the 
transverse channel approach smoothly the real axis as momentum increases. 



behaviour is also different from the poles found for gauge invariant operators. At weak coupling, 
these poles open up in branch cuts at fixed positions in the imaginary axis, while the holographic 
computation predicts that at infinite coupling the only singularities are poles that come closer to 
the real axis. Eventually a new peak appears in the spectral function, located close to ro = q. This 
peak persists at higher momentum, and can be interpreted as a quasiparticle excitation of charge 
density fluctuations. This shows a change of behaviour of the system as we increase the momen- 
tum, from diffusive to reactive. Transverse fluctuations are comparatively featureless, there are no 
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quasiparticle excitations appearing; this reflects the fact that there are no propagating modes in the 
zero temperature conformal theory. 

An application of our results is the calculation of the hydrodynamic time scale Th- Using our 
numerical results we find that in units of {2nT)^^ the minimal time scale is Th = 3.7 — 3.2 in a 
rank q = 0.3 — 0.48, rapidly growing for higher values of the momentum. For lower values it also 
grows, but this is due to the fact that the charge distribution is already quite uniform, so we can 
take this value as the onset of diffusion. In fact the dispersion relation starts to deviate from the 
hydrodynamic approximation for q ^ 0.45 which corresponds to Compton wavelengths of 1.2 fm, 
the size of a proton. As a model for sQGP in RHIC physics, let us pick T ~ 2Tc ^ 350 MeV. The 
hydrodynamic time scale is then Tr ~ 0.3 fm/c. This is a remakably short time. In fact it indicates 
that the hydrodynamic approximation is valid from very short times on. Notice that at RHIC the 
thermalization time is Ttherm ~ 0.6 fm/c and the hydrodynamic approximation is therefore valid 
already for t < I fm/c We expect that the values for shear or sound modes will be slightly 
different, but it is reassuring to find the right orders of magnitude even for !N=4. 
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